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Abstract 

The ground state of a one-dimensional spin-1/2 chain with periodical boundary condition in the 
Heisenberg XY model is investigated. We consider the spatial correlation and concurrence between 
any nearest-neighbor pair of spins under the conditions of different coupling strength, anisotropic 
parameter and magnitude of a transverse field. Quantum phase transitions due to the competition 
between coupling and alignment which cause the abrupt changes of the correlation and concurrence 
are observed. The transition are direct results of the level crossing. 
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I. INTRODUCTION 



The quantum Heisenberg XY model, which was first intensively studied by Lieb et al 
in the early 1960s []|], has aroused new interests in recent years. It was taken to be a 
potential candidate model in the practical apphcations in quantum computation [2|, y, |4| 
and quantum information j^, ^. The Heisenberg model is not only a theoretical model 
with many interesting physics in it, but also realized in recent years in laboratory in 
physical systems such^ quantum dots [2], nuclear spins [8] and electronic spins [9|]. Many 

had been devoted into the model to explore the irnportant 
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theoretical researchs 
quantum properties, especia 
quantum phase transition 



ly the correlation function 12(], the entanglement 13|] and 
[itI . [isl in one dimensional systems. The studies show 
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that the amount of the spatial correlation and entanglement between two spins in a spin 
chain (with open end) or a spin ring (with periodical boundary condition) can be modified 
by changing the magnitude of the external magnetic field and/or the interaction strength. 
And these two quantities shows singularity and obeys the scaling law in the vicinity of the 
quantum phase transition point of the system. 



The quantum phase transition, which can take place at T = 0, is an interesting subject 
and intensively studied in the last decades fl^ . 11^ . In this paper, we studied the quan- 
tum phase transitions of the XY model through the calculation of the correlation and the 



concurrence 
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21| between two nearest-neighbor spins. By changing the magnetic field or 



coupling coefficient, we find that both of the two quantities have finite jumps at some critical 



points. 

crossing 



he states of the system are carefully studied, the jumps are related to the level 
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24| of the ground state to the excited state. The change of the symmetry 



of the ground state [25|, l26|, |27|, l28|, l29|, l30| is the course of the quantum phase transition. 
The rest of this paper is organized as follows. In section [Tll we introduce the Hamiltonian 
for the Heisenberg XY model in an external field and describe the computation procedures 
for the correlation and the concurrence of the subsystem; Detailed results and discussions 
are given in section llllt The conclusion of our study is given in section IIVI 
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II. THE HAMILTONIAN AND THE THEORY 



The model studied here is a one dimensional spin-1/2 chain with periodical boundary 



condition (PBC) 
is: 
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33] under a transverse isotropic magnetic field. The Hamiltonian 

N N 



■t=l 



i=l 



where N is the number of spins, cr^ = cr^ ~ froi^ the PBC. Jx and Jy 

are the coupling strength between nearest-neighbor spins along the direction of x and 
respectively, is the amplitude of the magnetic field along z direction, a^., (jy and cr^ are 
the well-known Pauli matrix: 



1 

1 



, (yy 



-i 

1 



1 
-1 



(2) 



The parameters Jx and Jy may be represented by their average value J and anisotropic 
parameter 7 as: 

J, = (1 + 7) J, J, = (1-7) J, (3) 
and the Hamiltonian may be written as: 



N 



N 



(4) 



1=1 



1=1 



where a± = \{crx ± if^y)'- 
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0^ 









(5) 



The spatial correlation [3J, l35|, l36| between two spins, 1 and 2, is defined by 

Cl,2 



(1) 



(6) 
(7) 
(8) 



P12 = Tr'(p). 

Where a^^^a^"^^ can be represented explicitly in a matrix form in the product eigenspace of 
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(9) 



/l o\ 
0-120 
2-10 
\0 1/ 

And pi2 is the reduced density matrix of the spin 1 and 2 derived from the density matrix 
ps of the whole system. Tr'(- ■ ■ ) means to trace out all the spin degrees of freedom except 
spins 1 and 2. In this study, we are interested in the ground state of the whole system, 
which is the lowest energy eigenstate ipo of H, so ps = \ipo){i^o\- 



The concurrence of spins 1 and 2 is defined as [20|, l2l|] : 

Con = max{Ai — A2 — A3 — A4, 0}, (10) 

where Aj are the square roots of the eigenvalues of the product matrix puPu in decreasing 
order. Where pi2 is the reduced density matrix of the two spins subsystem, and pu is 
constructed as {ay ® o-y)p*{ay ® ay). Equation (fTOj) applies to both of the mixed state and 
pure state. The reduced density matrix may be write down implicitly as: 

/ Pl,l Pl,2 Pl,3 Pl,4 \ 

P2,l P2,2 P2,3 P2,4 

P3,l P3,2 P3,3 P3,4 

\ P4,l P4,2 P4,3 P4,4 / 



P12 



(11) 



then we get pi2 as the matrix below: 

/ 



P12 



P4,4 "~P3,4 ~P2,4 Pl,4 
~P3,4 P3,3 P2,3 ~Pl,3 
^P2,4 P2,3 P2,2 — Pl,2 

V P^,4 -Ph ~Ph Pl,l / 
If the bipartite quantum state pi2 is pure, it can always be written as: 

P = 

1^) = a|00) + 6|01) +c|10) + c/|ll), 
then equation flTUl) could be simplified to 

Con{\il))) = 2\ad-bc\. 



(12) 



(13) 
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TABLE I: The discontinuities in Fig. [T] 



N 


J (discontinuities) 


4 


0.649 


1.569 


/ 


/ 


/ 


6 


0.650 


0.888 


2.426 


/ 


/ 


8 


0.650 


0.767 


1.148 


3.268 


/ 


10 


0.650 


0.721 


0.908 


1.415 


4.104 



TABLE IL The discontinuities of the case of = 6 in Fig. [2] and H 



Bz 


J (discontinuities) 


0.65=1.30 X 0.5 


0.325=0.650 X 0.5 


0.444=0.888 x 0.5 


1.213=2.426 X 0.5 


1.30=1.30 X 1.0 


0.650=0.650 X 1.0 


0.888=0.888 x 1.0 


2.426=2.426 x 1.0 


2.60=1.30 X 2.0 


1.300=0.650 X 2.0 


1.776=0.888 X 2.0 


4.852=2.426 x 2.0 



III. RESULTS AND DISCUSSIONS 

First, we consider the isotropic case, 7 = 0, in an external field = 1.30. Figure [1] 
gives the correlation functions as function of coupling strength J for different sizes of the 
system, (A^ = 4, 6, 8, 10). From the figures, we see that with the increase of J, C12 displays 
a "ladder" behavior. We also noted that for the case of spins in the system, there are 
N/2 discontinuities, the position of the jump can be found in table HI 

The effect of the magnitude of magnetic field on the phase transition was shown in 
figure [21, in which we plot the case of a 6 spins system with three different magnetic fields, 
Bz = 0.65, 1.30 and 2.6. The corresponding transition points of J are listed in table [III 
By simple calculation and comparison, we found that the transition values of J for each 
transition is approximately proportional to the magnitude of the external field. 

Figure [3] plots the variation of correlation C12 with the applied field B^ for fixed coupling 
coefficient J = 1.0 and different values of anisotropic parameter 7. The system sizes 
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are = 6. When I7I < 0.8, there are 3 jumps for a 6 spins system as in the case of 
isotropic case. However, with the increase of anisotropy 7, the magnitude of the jumps 
becomes smaller and the transition points moves slightly to the direction = 0, till finally 
disappears in the completely anisotropic case 7 = 1. Since the 7 = 1 corresponding to the 
Ising case, where the internal number of freedom is 1, the phase transition is different to 
the XY case where the internal number of freedom is 2. Our calculation shows that the 
cross over from n = 2 to = 1 is continuous. It is also verified that there is no quantum 
phase transition in the quantum Ising case. 

Now we turn to the concurrence between two spins, Coni 2- The following results were 
obtained for a system of = 6. Figures H] plot the variation of concurrence Coni^2 with 
the coupling strength J in the isotropic case (7 = 0), the applied magnetic field fixed to 
three values, = 0.65, B^ = 1.30 and B^ = 2.60. From the figures we see the same 
phase transition behavior as that of the correlation functions. Jumps are found at exactly 
the same places, which can also be checked in table [Tll as for the correlation functions. 
Figures [5] gave the variation of Coni^2 with the magnitude of external field Bz for 8 different 
anisotropy 7. When I7I < 0.5, the jumps in concurrence are clearly observed while for larger 
|7|'s the magnitude of the jumps becomes smaller and smaller and turn to the smooth curve 
in the case of 7 = 1. 

It is interesting to note that the physics behind the above phase transitions turns out to 
be very simple and clear, it is the result of level crossing of the ground state. To clarify this 
point, we turn to a detailed analysis of the ground states of the system under different sets 
of parameters. To be specific, we use |1) and |0) to represent the "spin down" and "spin 
up" state for each spin in the eigenspace of {cr^,^^} respectively. The state of the whole 
system of spins ip can be written as a linear combination of 2^ states and 2^ normalized 
coefficients: 

^ = ao|00---00) + ai|00---01) + --- + a„|ll---ll), (14) 

n 

n = 2^-l, ^|a,|2 = l. 

where |00---00) stands for the direct product state |0)i ® ■•■|0)7v We diagonalized the 
Hamiltonian and calculated the eigen function and the corresponding eigen energy. We take 
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the coupling strength J = 1 and anisotropy 7 = for simphcity. Figure [6] shows the variation 
of several states energy (Someone of them will be the ground state in some interval of Bz) 
with the magnitude of external field B^. It is clear from the figure that level crossing occurs 
at the positions where abrupt changes of the correlation function and the concurrence were 
observed. Concretely, when < -B^ < 0.525, the ground state (which has been normalized) 
of the whole system is 

1^1 =~ 0.11785|000111) + 0.23570|001011) - 0.23570|001101) + 0.11785|001110) 

- 0.23570|010011) + 0.35355|010101) - 0.23570|010110) - 0.23570|011001) 

+ 0.23570|011010) - 0.11785|011100) + 0.11785|100011) - 0.23570|100101) (15) 
+ 0.23570|100110) + 0.23570|101001) - 0.35355|101010) + 0.23570|101100) 

- 0.11785|110001) + 0.23570|110010) - 0.23570|110100) + 0.11785|111000). 

Which is basically the combination of the states with antiferromagnetic order, that is, the 
states with half of the spins "up" and half of the spins "down". When 0.526 < B^ < 1.464, 
the ground state changed to 

^0= + 0.16667|000011) - 0.28868|001001) + 0.16667|001010) + 0.33333|001001) 

- 0.28868|001010) + 0.16667|001100) - 0.28868|010001) + 0.33333|010010) 

- 0.28868|010100) + 0.16667|011000) + 0.16667|100001) - 0.28868|100010) 

- 0.28868|101000) + 0.16667|110000). 

This is the state that 2 spins "up" and 4 spins "down". In the course of increasing the 
applied field B^, the energy level corresponding to the state (fT5|) increases and the energy 
level corresponding to the state f|T6|) decreases. At the critical field, B^ = 0.525, the two 
levels meet and then crossed, then the ground state of the spin chain is a combination of the 
two level states. Beyond the critical field the ground state becomes the state ( |T6l) . The cross 
of the ground state changed the symmetry of the ground state abruptly so that a quantum 
phase transition occurs. At about Bz = 1.465, another level cross take place and over that 
point, the ground state of the system becomes, 

^3 = -0.40825|000001) + 0.40825|000010) - 0.40825|000100) + 0.40825|001000) 
-0.40825|010000) + 0.40825|100000). 

This is a combination of states with 1 spins "up" and 5 spins "down". And finally, when 
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Bz > 2.001, the ground state is simply as 

jp^ = 1 000000), 

which is the state with ferromagnetic order. We see from the above analysis that the 
quantum phase transition is a result of the level crossing of the ground state, and the level 
crossing is simply the result of the competition between the coupling strength J, which 
favors antiferromagnetic order, and the applied field, which favors the ferromagnetic order. 
By increasing the applied field while keeping the coupling strength constant, the ground state 
changed from antiferromagnetic order to the ferromagnetic order in several steps, turns a 
spin from "up" to "down" at each step. So we observed that V^o ^ V'o V^o "^o- This 
also explains the number of jumps in the correlation function and concurrence, which must 
be half of the whole number of spins. 

IV. CONCLUSION 

In this paper, we studied the spatial correlation Ci,2 and the degree of entanglement 
Coni 2 between two nearest-neighbor spins along a one-dimension spin ring in Heisenberg 
XY model. We found that in the isotropic case with fixed external field B^, quantum phase 
transitions of both Ci 2 and Coni^2 were observed with the increment of coupling strength J, 
the phase transitions will smoothed out as the anisotropy 7 increased and disappears at the 
extremum value of 7 = 1. The transition is a direct result of the level crossing of the ground 
state of the whole system, in which the symmetry of the ground state changed abruptly at 
the crossing point. 
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Figure Captions 



FIG.d) The nearest-neighbor correlation Ci,2 in different size of the spin ring as function 
of couphng strength J, parameters are 7 = and Bz = 1.30. 

FIG. |2) The nearest-neighbor correlation Ci^2 under three different cases of external 
field: = 0.65, B^ = 1.30 and B^ = 2.60, the parameter 7 = 0. 

FIG.m The dependence of nearest-neighbor correlation Ci 2 on the applied field B^ with 
different anisotropic parameters 7. All the results are calculated with J = 1.0. 

FIG. m The comparison of Concurrence for two spins in the isotropic case, 7 = 0, in 
three different applied field: Bz = 0.65, B^ = 1.30 and B^ = 2.60. 

FIG. [5l The dependence of Concurrence Concurrence for two spins on the applied field 
Bz with different anisotropic parameters 7. All the results are calculated with J = 1.0. 

FIG. [6l The energy levels as function of Bz for given coupling strength, the system size 
is = 6, some of the levels which do not take part in the level crossing were omitted for 
clarity. The other parameters are: 7 = 0.0, J = 1.0. 
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